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Abstract 

Electrodynamic phenomena related to vortices in superconductors have 
been studied since their prediction by Abrikosov, and seem to hold no fun- 
damental mysteries. However, most of the effects are treated separately, 
with no guiding principle. We demonstrate that the relativistic vortex 
worldsheet in spacetime is the object that naturally conveys all electric 
and magnetic information, for which we obtain simple and concise equa- 
tions. Breaking Lorentz invariance leads to down-to-earth Abrikosov vor- 
tices, and special limits of these equations include for instance dynamic 
Meissner screening and the AC Josephson relation. On a deeper level, we 
explore the electrodynamics of two-form sources in the absence of elec- 
tric monopoles, in which the electromagnetic field strength itself acquires 
the characteristics of a gauge field. This novel framework leaves room for 
unexpected surprises. 



1 Introduction. 

The study of the matter formed from Abrikosov vortices in type-II supcrconduc- 
tors[l| constitutes a vast and mature research subject. This subject is crucial for 
the technological applications of superconductivity^ but it has also proven to be 
a fertile source for fundamental condensed matter physics research. The clastic 
and hydrodynamical properties of matter formed from vortices can be very easily 
tuned by external means and it has demonstrated to be an exceedingly fertile 
model system to study generic questions regarding crystallization, the effects of 
background quenched disorder and so forth0, 0] ■ Especially after the discovery 
of the cuprate high-T c superconductors it became also possible to study the 
fluids formed from vortices. Because of the strongly two-dimensional nature 
of the superconductivity in the cuprates, the Abrikosov vortex lattice becomes 
particularly soft and it melts easily due to thermal motions at temperatures 
that are much below the mean field 7J C 2-line 

Many phenomena in this field are of a dynamical nature, associated with the 
fact that vortices are in motion. This includes the vortex flow, the magnetic 
field penetration and the flux creep, but also the large Nernst effect of the 
vortex fluid and, perhaps most spectacularly, the use of cuprate vortices as 
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source of terahertz radiation [|| 0|. This vortex dynamics is analogous to the 
magnetohydrodynamics of electrically charged plasmas in the sense that the 
forces exerted on vortices are exclusively of electromagnetic origin, while in turn 
the vortex matter backreacts on the electromagnetic fields. The phenomena 
that arise are rather thoroughly understood departing from the AC-and DC- 
Joscphson relations as well as the Maxwell equations as the force equations in 
this "vortex magneto-hydrodynamics" . 

Although the computations explaining these phenomena are correct, they are 
of a rather improvised ad hoc nature, at least compared to the Landau-Lifshitz 
styleQ of deriving the usual magnetohydrodynamics from first principles. In 
so far as the forces are concerned, the latter eventually departs from the micro- 
scopic action describing Maxwell electrodynamics, 

S = J d 4 X /^Maxwell, (1) 
-^Maxwell ~^F^ V F^ U -|- A^J^ L: (2) 

where the conserved currents J M parametrize the world lines of the charged 
particle forming the plasma, sourcing via the electromagnetic gauge fields 
the electromagnetic field strength F M „ = d^A v — dyA^. Varying the action with 
respect to A v one obtains the Euler-Lagrange equations of motion, 

dfiFfiv = Jvi (3) 

the Maxwell equations in relativistic shorthand notation. What are the corre- 
sponding equations, dealing with the Abrikosov vortices sourcing the electro- 
magnetic fields? In a related but yet different context Q we stumbled on this 
question, finding out that the answer is apparently not available in the literature. 
To our impression this is rooted in the fact that one needs a piece of mathe- 
matical technology that is unfamiliar to condensed matter physicists, while it 
is well known in differential geometry [l^ and string theory [Tl|. We refer to the 
formalism of two-form gauge theory. This amounts to a generalization of the 
familiar "one-form" Maxwell (and Yang-Mills) gauge theories characterized by 
gauge fields (vector potentials) A^ that carry one label. Instead one has gauge 
fields b^ u , representing antisymmetric tensors, that have gauge transformations 
in the form of the addition of the gradient of a vector field (that therefore looks 
like a Maxwell field strength = <9 M a„ — <9„a M even though the meaning is 
completely different): — > b^ + f^. How can such a structure arise? 

The key ingredient is that vortices as topological excitations of the supercon- 
ductor are lines in 3-dimensional space when they are static. Invoking time as 
one should using the action principle, these lines spread out in space time as sur- 
faces or world sheets, in the same way as point particles correspond with world 
lines. In fact, in string theory, Abrikosov vortices are known as "Nielsen-Olesen 
strings" [l2j , regarded as rather primitive as compared to the fundamental "crit- 
ical" strings. We can now rely on Schwinger's principle [l3| that the sources are 
the principal objects, emitting and absorbing the force-mediating gauge fields. 
Dealing with point particles as the sources, their infinitesimal world line inter- 
vals are parametrized in terms of a charge density p and current J which can be 
combined in the spacetime covariant current = (cp, J) of Eq.'s @ and ([3]), 
that satisfies a continuity equation, 

0" J M = dtp + V ■ J = 0. (4) 
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The continuity equation is also a constraint that removes the longitudinal com- 
ponent of J M . The fields that couple to J M are therefore gauge fields, since 
the longitudinal component is not sourced and as a consequence not physical. 
In this context of electromagnetic currents, the field just corresponds with 
the photon gauge field. 

In the case of the world sheets formed by the vortices in 3+1-dimcnsional 
spacetime one has to depart instead from infinitesimal world sheet areas and 
these have to be parametrized in antisymmetric tensors J^ v . In analogy with the 
particle currents, these "two-form currents" will source two-form gauge fields, 
subjected to the transversality condition. The vortex action has therefore to 
contain a term b^ v J^ v where b^ v is a two- form gauge field. As we will show, 
for the vortices in 3+1 dimensions the pieces fit in a remarkably elegant way: 
instead of the gauge potential A^, now the physical electric and magnetic fields 
E and B as collected in the Maxwell field strength F^, couple minimally as 
two-form gauge fields to the vortex world sheet current. As we will show, this 
gauge coupling involves actually the (Hodge) dual of the electromagnetic field 
strength, so that the coupling term is F^e^" J^ A , where is the fully 

antisymmetric Lcvi-Civita symbol. Remarkably, in so far this "Schwingcr side" 
is concerned, the one-form Maxwell theory is lifted by an emergent gauge invari- 
ance to a two-form gauge theory. Given that inside the type-II superconductor 
only vortices are present as sources for the electromagnetic fields, and since these 
carry only magnetic fluxes and no electrical monopole charges, the Maxwell field 
strengths of the fundamental vacuum lose their physical meaning and turn into 
gauge fields. However, there is a caveat: besides the vortices the electromag- 
netic fields also experience the Meissner- and Thomas-Fermi screening. These 
restore the physical reality of the Maxwell field strengths, explicitly breaking 
the two-form gauge invariance. Only in the extreme type-II limit, where both 
the Meissner and Thomas-Fermi screening lengths diverge, the emergent two- 
form gauge invariance becomes literal. Otherwise the two-form gauge theory is 
forced in a "fixed frame" corresponding with the generalized Lorenz gauge fix 
of the fully gauge invariant theory. This peculiar gauge theoretical structure is 
the key to the special nature of vortex electrodynamics. 

On the practical side, with this mathematical technology at hand it is 
straightforward to derive the fundamental action and equations of motion gov- 
erning the electrodynamics of vortices in type-II superconductors: Eqs. (|16|) and 
(|5T|l , (|64[) & ([S5]l for the simpler relativistic- and the realistic non-relativistic vor- 
tices, respectively. These originate from an action Eq. ((29)) resp. (|56|) . Merely 
for reasons of convenience we will develop much of the theory in the relativistic 
limit where it is assumed that the superfluid phase velocity is coincident with 
the speed of light. As in electrodynamics, this greatly simplifies the equations 
and it is easier to follow the argument. Of course, in real superconductors this 
phase velocity is a small fraction of the light velocity and it is straightforward 
to break the Lorentz invariance to obtain the equations presented in Section [6] 
which are of relevance to earthly superconductors. 

We claim that these equations of motion represent a fundamental and com- 
plete description of the electrodynamics of the Abrikosov vortices in super- 
conductors. We will demonstrate that all known physical properties associated 
with this electrodynamics are economically encoded in this formalism: Thomas- 
Fermi screening and dynamical Meissner screening, the electric fields induced 
by vortex motion and the Ncrnst effect, the AC Joscpson effect and radiation 
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due to vortex motion. These all arise as special limits of our general equations. 
Admittedly we have not achieved extracting new vortex physics but we believe 
that there is a potential for surprises. We invite the expert readership to have a 
closer look to find out whether they can profit from this particularly convenient 
mathematical formulation of the problem. 

The paper is organized as follows: the fundamental definition of a vortex 
is reviewed in Section [21 and it is shown how vortices are described as world 
lines and world sheets. In Section[3]we show how, in a superconductor, the vor- 
tex world sheet dynamics follows directly from the time-dependent Ginzburg- 
Landau equations by a single mathematical derivative operation. This is the 
easiest way to actually obtain our central results. The next two sections have a 
theoretical emphasis: in Section 0] we motivate the general mathematical struc- 
ture of the electrodynamics sourced by two-form fields. We demonstrate here 
that the field strength F M „ becomes itself a gauge field. A rigorous deriva- 
tion of the action governing the electrodynamics of vortex sources using duality 
techniques is presented in Section [5] This comprises the dual gauge field that 
mediates interactions between vortices. Finally, by breaking the Lorentz in- 
variance we derive the equations of motion of vortex electrodynamics from this 
action in their full extent in Section [51 and we derive the physical phenomena 
mentioned in the previous paragraph. Finally, in the Appendix this structure is 
formulated in the general mathematical language of differential forms, allowing 
us to extend the results to any spatial dimension larger than 2. 



2 Vortex fields. 

Let us consider vortices in type-II superconductors. These are of course the 
familiar flux lines where the magnetic field penetrates through the supercon- 
ducting sample, but the more profound statement is that a vortex is a region 
where the superconducting phase is singular. Let us repeat precisely what is 
meant by this for the case of a 2+1-dimcnsional superconductor. This argument 
applies just as well to a supcrfluid (uncharged superconductor). More details 
are found in Kleinert's textbook [3, ch. II. 1]. 

In the Ginzburg-Landau description, the parent state of a superconductor 
obeys a t/(l) symmetry, as the phase ip of the complex order parameter \I/ = 
|^/|e lv is only defined up to global transformations \& — > Vfe 1 ". However, across 
the phase transition, the amplitude obtains a vacuum expectation value, and 
one of all possible phases is singled out. This is the well known phenomenon of 
spontaneous symmetry breaking. Hence, the superconductor is a state where 



the phase has a preferred value over large length scales (Fig. 1(a)). It is now 
possible that the system is disturbed so that a vortex is formed. In two spatial 
dimensions, a vortex is pointlikc (a 'vortex pancake'). Tracking the deviation 
of the phase from its preferred value along a contour C around the vortex origin 
reveals that the phase changes by 2irN, where N is the winding number (Fig. 



1(b)). The mathematical statement is, 

d^ip = 2nN. (5) 



1 ln a superconductor this is lifted to the local U(l) gauge symmetry of electromagnctism, 
but that does not matter here as is well known. As wc will see explicitly, the gauge symmetry 
can be included in a next step. 
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(a) P hase coherence . S pont a- 
neous symmetry breaking singles 
out a preferred phase value. 




(b) Vortex. A single vortex of winding num- 
ber 1. Tracking the phase change along the 
contour C (red) adds up to 2iz. Also shown in 
hatched red is the surface S enclosed by C. 



Figure 1: Phase coherence and its distortion by a vortex. The hallmark of the 

superconductor is the long-range correlations in the order parameter phase. The phase is a com- 
pact variable, and as such permits vortex configurations, that cause the phase to change along its 
perimeter in multiples of 1-k only. 



We can now formally invoke Stokes' theorem, by converting the contour integral 
into a surface integral of the curl of the integrand, over the surface S enclosed 
byC, 

/ dS x e x ^ d v d^ = 2-kN. (6) 
Js 

This is allowed rigorously only when ip is non-singular in all of S, but we turn 
necessity into virtue by noticing that the above equation is satisfied when, 

e Xvfl d v d^ip{x) = 2TrN6^\x). (7) 

Here 5 X (x) is a two-dimensional delta-function in the plane orthogonal to A. 
It is only non-zero at the origin. This corresponds to the fact that for a vortex 
<p is singular at the origin. For the same reason, the derivatives on the left-hand 
side do not commute and the expression is not manifestly zero, as one would 
expect for the contraction with the fully antisymmetric tensor e\ v ^ . 
We define, 

Jj (x) = e\ ull d u d^<p(x), (8) 
as the vortex current. Its temporal component, 

J t v (x) = (d x d y - dyd x )tp{x) = <5 (2) (x), (9) 

is the density of vorticity, just as p is the charge density related to the electro- 
magnetic four-current J M = (cp, J). Similarly, the spatial part J ; v represents the 
current or the flow of the vortex. Together, the density and the current satisfy 
a continuity equation, 

d x jy = a t j t v + d t jy = o. (io) 

The density of vorticity can only decrease (increase) when the vortex flows away 
from (towards) the current location. In terms of Eq. ([7]) , the continuity equation 
is a conservation law or intcgrability condition, signalling that vortex world lines 
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(a) Vortex world line. Segment of a vortex world 
line (red) in 2+1 dimensions. Shown in blue is the 
line element J ^ and it temporal and spatial com- 
ponents Jt and J. For an electric monopole, these 
correspond to the charge and current densities p 
and J. 




(b) Vortex world sheet. Segment of a 
vortex line tracing out a world sheet in 
time. The vortex lies in the x-y-planc, 
the third spatial dimension cannot be 
drawn. Shown in blue is the world sheet 
surface element J K \. 



Figure 2: Vortex WOrld line and WOrld Sheet A graphical picture of dynamical system 
is provided by tracing the motion of particles or line objects through time. The traces arc called 
world lines or world sheets. In field theory, they are represented by one- or two-form fields, which 
represent the line or surface elements of the world line or sheet. If the world lines or sheets can not 
end in the middle of a system, the one- or two-form fields are conserved, as cLJ/i — or d K J K \ — 0. 



cannot begin or end in the middle of a superconductor. It is revealing to regard 



as representing the world line of a vortex itself (Fig. 2(a) ). At each point in 
spacetime x, J\(x) is the line element of the world line; its temporal component 
reflects the density (in multiples of 2tt) and the spatial components indicate in 
which direction the vortex is moving. 

With the physical meaning of the mathematical definition Eq. ([8j) being 
explicit, we are now ready to generalize to 3+1 dimensions. In three spatial 



dimensions the topological defects are vortex lines. In spacetime such a line 
traces out a world sheet. In direct analogy with J^(x) being understood as the 
line element of a vortex world line, we now define, 

J K v A (x) = e KXvii d v d^{x) = 2nN8%(x), (11) 

as the surface element of the vortex wor Id sheet (Fig. [2(b)) . In four space-time 
dimensions, the Levi-Civita symbol € K x Ufjt has four indices. The vortex current is 
antisymmetric in its indices «A, which should be the case since those indices can 
never be parallel when we are to define a surface element. The delta- function is 
in the plane perpendicular to the world sheet. 

What is the physical meaning of the six independent components of JY\t 
Analogous to the situation in 2+1 dimensions, the temporal components J 4 Y 
represents the density of vorticity of a vortex line along the /-direction, while 
J^i is the flow of vorticity in the fc-direction of a vortex line along the /-direction. 
Since vorticity is topologically conserved, there cannot be any vortex flow in the 
direction in which it is pointing, hence the antisymmetry in k and I. For each 
spatial direction / there is a continuity equation relating the decrease (increase) 
in vorticity to its flow away (towards) the current location. This may be sum- 
marized as d K J^ x = 0. This generalizes the continuity equation to the flow of 
strings; notice that this involves three independent conditions instead of the sin- 
gle condition one finds for particles. To gain some intuition, imagine a straight 
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vortex line in the z-direction, the density of which is given by J^ z . It can move in 
either the x- or the y-direction, denoted by J^ z and Jy Z respectively. Together 
they satisfy the continuity equation, 

d t jZ + d x Jl + d y jv z = o. (12) 

In the following sections it will become clear that treating vortex dynamics 
in terms of its world sheet dynamics is simple and straightforward, allowing for 
immediate derivations of many standard relations in superconductivity. 

3 The vortex world sheet in relativistic super- 
conductors. 

We will now show how the vortex world sheet appears from the Ginzburg- 
Landau equations. In section 21 we shall derive the more generic coupling of a 
vortex current to electromagnetic fields. 

Since we are pursuing a relativistic treatment, we adopt an action rather than 
a Hamiltonian approach, where time coordinates and derivatives are treated on 
equal footing with their spatial counterparts. At any time, the equations familiar 
from condensed matter can be obtained from the spatial part of the Lagrangian. 

Before we write down the partition function let us stress that it may be less 
familiar to researchers in the field of superconductivity, since it will be fully rela- 
tivistic. In particular it will have a squared time-derivative, whereas most works 
start with a single time-derivative term. The latter applies to systems which 
are diffusion-limited. Of course, in actual superconductors vortices are accom- 
panied by such diffusion processes. However, the relativistic action is necessary 
to derive the vortex worldsheet. Futhermore processes such as Thomas-Fermi 
screening are in fact ballistic. Finally the validity of this relativistic approach is 
verified by the results of Section [51 If one wishes to consider diffusion processes, 
an appropriate term can be added to the Lagrangian at will. 

The partition function associated with the relativistic Ginzburg-Landau ac- 
tion deep in the superconducting state is, 

Z = j V<pVA r F(A^)<il h (13) 

(14) 

Here F^ u = d^A v — d u A^ is the electromagnetic field strength; denotes 
an appropriate gauge fixing condition; <p is the superconducting phase related 
to the order parameter = ^/p^c lip ; p s is the superfiuid density; m* and e* are 
the mass and charge of a Cooper pair; and most importantly, one must take 
great care to differentiate between the two velocities in the problem, namely 
the velocity of light c pertaining to the photon field A^, and the phase velocity 
in the superconductor c p h. Therefore we have defined = (do, V), do = ^.dt 
and 9P h = (<9jf\V), d? h = ±d t . Furthermore A^ = (-\V,A) and A^ h = 
(—■^V, A). The last form is dictated by gauge invariance of the second term 
in Eq. (|14p. The metric is rf" v = diag(— 1, 1,1,1), but we shall use only lower 
indices for notational simplicity. 
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We shall for the moment proceed in the relativistic limit where c p h = c, for 
simplicity. The equations of motion then follow from variation with respect to 
A u , 

dC dC 1 h 2 e* e* 



^{^KXpud^dpA,, — CjiXppd^d^dnAp) + (. K \ pv d p A v — —e K \p V d p d v ip — —J^\- 



Now we act with t K \ pv d p on this equation, which leads to, 

t \pvdpd v ip = 

(16) 

Here we have dchncd the London penetration depth A = y ^ "*" 2 p ~' the second 
term vanishes because the antisymmetric contraction of two derivatives; and 
on the right-hand side we recognize from Eq. ([TTjl the definition of the vortex 
current J^ x . Let us consider the special case k = t, and use the definition of 
the magnetic field B\ = ei rn d r A n , 

- X^B, + Bi = A jV = t»! 2) (4 (17) 
f* e* e* 

Here we have used relation Eq. (fTTj). This is precisely the textbook equation 
for the Meissner screening of a vortex source of strength TV, with flux quantum 



$o = 27rft/e* [la, eq.(5.10)]. But instead of ad hoc inserting the delta-function 
source, wc actually derived it from the singular phase field. The only difference is 
that here also the dynamics are taken into account via the double time derivative 
contained in oz. The true power of the vortex world sheet shows itself when 
considering the electric field E = — V^4o — dtA. and the spatial components 
of the vortex field. This will be further elaborated on in Section [6] But let 
us first analyze how two-form sources couple to clcctromagnctism in general, 
followed by a more general derivation of the above relations invoking a duality 
mapping, by which wc can treat the vortex fields in the action itself, rather 
than only in the equations of motion. This can be regarded as revealing the 
more fundamental structure of the problem. The reader who is less interested 
in these theoretical matters may skip ahead directly to Section [5] 



4 Electrodynamics of two-form sources. 

We will formulate here the generalization of the standard Maxwell action and 
equations of motion when the sources are not monopolcs with charge density p 
and current J m , collected in a vector field J M = (cp, J m ), but instead (vortex) 
lines with line densities J t i and line currents J^i (which denote the current 
in direction k of a line that extends in direction I), collected in a two-form 
field J K \ — (Jti, Jki)- Let us first recall the established knowledge for ordinary 
electromagnetism, in terms suited for this generalization. For clarity reasons we 
again use a shorthand notation where we are intentionally sloppy with contra- 
and covariant indices, leaving out dimensionful parameters in order to maximally 
expose the principles. In the next section we will present the final results that 
are accurate in this regard. 

Let us start considering a set of electrical monopole sources collected in a 
source field J M as in the above, satisfying a continuity equation/conservation 
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law c^Jp = 0. These sources interact via the exchange of gauge particles, as 
gauge fields that couple locally to the source fields, by an interaction term 
in the Lagrangian of the form A^J^. Because of current conservation, any 
transformation of the gauge field A^ — > A^ + d^e, where e is any smooth scalar 
field, will leave the coupling term invariant. Indeed, 

AfxJfj, -> A^ + (<VK* = A^Jp - ed^ = A^. (18) 

Here we performed a partial integration in the second step. The field strength 
F/u/ = dp.Au — d u A^ is also invariant under the same gauge transformation. An 
immediate consequence of this definition are the Bianchi identities or homoge- 
neous Maxwell equations, 

CaPuvdpFnv = taPnudpdpAv = 0, (19) 

because the derivatives commute. These equations comprise V • B = and 

V x E = — dfB. This suggests a Lagrangian of gauge invariant terms, 

^Maxwell = —-^FnvFfa, + Apjp, (20) 

accompanied by the Euler-Lagrange equations of motion obtained by variation 
with respect to A v , 

dpF^u = -J v . (21) 

These are the inhomogeneous Maxwell equations comprising V • E = p and 

V x B — <9 t E = J. In a superconductor, one must also add a Meissner term, 
which in the unitary gauge fix turns into a mass term for the gauge field A^, 

^Maxwell + Meissner — ~tF^ i i/F^ i1 i —A^A^ -\- A^J^^ (^2) 

In this form, the Meissner term breaks the gauge invariance of the Lagrangian. 
This corresponds to releasing the longitudinal degrees of freedom of the photon 
field. A gauge equivalent perspective is that this degree of freedom represents 
the phase mode of the superconducting condensate. The equation of motion is 
modified to, 

d„F^ -A v = —J v . (23) 

Let us now repeat this procedure for antisymmetric two-form sources J K \ = 
(Jti, Jki)- These must obey the continuity equations/conservation laws d K J K \ = 
0, reflecting that the density of the source can only increase (decrease) when it 
flows into (out of) the region under consideration, and that vortex lines cannot 
end within in the system (no monopolcs). Consider now that these sources 
interact by exchanging two-form gauge fields, that we will tentatively denote by 
G K \. Then these gauge fields couple locally to the sources as G K \J K \. These 
fields have to transform under gauge transformations as, 

G K x^G KX + ^(d K E X -d x e K ), (24) 

where £\ is any smooth vector field, in order to leave the coupling term invariant 
as required by the current conservation. Indeed, 

G K \J K \ —> G K \J K \ + (d K e\)J K \ = G K \J K \ — £\d K J K \ = G K \J K \. (25) 
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Here we have used the antisymmetry of J K \ in the first step, and partial integra- 
tion in the second. The field strength H^ K \ = d^G^x] = dnG K \+d\G llK ,+d K Gx IJl 
is also invariant under these gauge transformations. An immediate consequence 
of this definition is the Bianchi identity, 

^v^n\dyH^ K \ = d^d ll G K x] = 0, (26) 

because the derivatives commute. With these definitions, we can write down a 
gauge invariant Lagrangian, 

C = -—H^ kX + G k xJ k x- (27) 

Note that this Lagrangian is in terms of the dynamic variables G K x, which we 
will see later is the dual of the electromagnetic field strength F^ u . In other words, 
this Lagrangian is in terms of the electric and magnetic fields themselves, rather 
than the gauge potential A^. The equations of motion follow after variation with 
respect to G K >, 

d^H^x = —JkX- (28) 

Now, in a gauge-invariance breaking medium such as a superconductor, one 
must add a Higgs or Meissner term to the Lagrangian as, 

£ = --^ h Ik\ - 4 G kA + G K xJ K \- (29) 

Up to now we have just reviewed the standard derivation of non-compact U(l) 
two-form gauge theory. Let us now specialize to the case of a vortex line in a 
superconductor. For such an Abrikosov vortex, we know that the density J^f is 
proportional to the magnetic field, and that the magnetic field is parallel to the 
spatial orientation of the vortex line. In fact, when the magnetic field intensity 
coincides with the lower critical field H&, the dimensionful vortex density may 
be denoted as before, combining Eqs. (fTTj) and (fTTl) . 

■# = <Mf°(r), (30) 

where $o is the flux quantum \. Because of these considerations, the vortex line 
density should couple to the magnetic field Bi. The definition of the Maxwell 
field strength is, 

Ftn = F n F mn — ^rnnl^li (^-0 

If we contract the last definition with ~^2 mn t-tbmm orLC finds Bi — eti mn F mn = 
Gu- Here we introduce the Hodge dual of the Maxwell field strength G a/ 3 = 
\^ap^vF^ v . Then the coupling of the vortex line density to the magnetic 
field Bi is written as GuJ^ and generalizes to G K xJ%x- Therefore, the general 
two-form gauge field in Eq. (|27p is now identified as the dual Maxwell field 
strength G K x- 

This has one immediate astonishing consequence: the Maxwell field strength 
F^u itself has now become a gauge field ! The gauge transformations Eq. (|24[) 
correspond to, 

Fp.v — > F^v + e^ VK xd K ex- (32) 
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How does it come about that these all too physical F^s have suddenly turned 
into gauge variant quantities? The reason is simple although perhaps defeat- 
ing the physical intuition: in normal matter we always have electric monopolc 
sources J v with the associated equations of motion d^F^ = — J u . In the ab- 
sence of any such sources, these equations reduce to d^F^ = 0. Together with 
the inhomogeneous Maxwell equations Cap^dpF^ = 0, these imply that the 
field strength cannot be measured at all. It amounts to the Schwinger wisdom 
that fields which cannot be sourced do not have physical reality. The formal ex- 
pression of this fact is that the field strength becomes pure gauge in the absence 
of monopole sources. 

Another insight is obtained by taking a closer look at the gauge transforma- 
tions Eq. (f3"2j) . For the Bianchi identities in Eq. (|19p these imply, 

tapfivdpF^ — > eapuvdpF^v + eaPtiudpefiunxd^ex 

= tctPiivdpFnv + ( d ad x - d 2 5 a \)e\. (33) 

In other words, the Bianchi identities are not invariant under these gauge trans- 
formations! This makes sense: these identities are a direct result of expressing 
the field strength in terms of a gauge potential A v , which of itself has three de- 
grees of freedom (four minus one gauge freedom). The Bianchi identities serve 
to restrict the six degrees of freedom contained in F^ v to the proper number of 
three. In the derivation of the two- form action Eq. (|?7|) . we have not assumed 
anything about the origin of the two-form field. Next to three physical degrees 
of freedom, there are three gauge degrees of freedom. Therefore the constraints 
e a piivdpF^ v = are not strictly enforced, but can always be obtained by a 
suitable gauge transformation. 

We never observe the gauge character of the fields F^ y themselves because 
the only two-form sources to which this action applies that we know of are 
Abrikosov vortices in a superconductor. The superconducting matter causes a 
finite penetration depth A for the fields, which is reflected by the addition of 
a Meissner term to the Lagrangian. The gauge-invariant form of this term is 
known to be, 

1 S d 2 - 8 8 

in the same way as one can formally write the Meissner term in Eq. (f212j) as 
F^v-rpF^y. However, since the longitudinal components of G K \ are not sourced 
by the conserved Abrikosov vortices, we are naturally led to the Lorenz gauge 
condition 8 K G K \ = 0, and the gauge freedom has been removed. With this 
gauge condition Eq. (f34|) reduces to G K \G K \, that appears in Eq. ((29|) . In 
other words, the superconducting medium forces us to the fixed frame action 
Eq. (HID. 

Wc end up with the action Eq. (|29l) . and we now put in dimensionful pa- 
rameters. Please note that this action is equivalent to the regular action as 
(|14p . but with the important difference that here we work with the dual field 
strength G K \ as the dynamic variable instead of the gauge potential .The 
equations of motion ( "Maxwell equations for relativistic vortices" ) are now ob- 
tained straightforwardly by varying with respect to G K \ as, 

A 2 (9 2 G A1!y — dp,d K G KU + dvd K G K ^) — G^ v = — ~J^ V - (35) 
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This is to be compared with Eq. ([25]) and Eq. (fTT)) . The second and third 
term can be set to zero by a gauge transformation Eq. (|32p or alternatively by 
invoking the Bianchi identities Eq. (|T5| . The meaning of these equations is that 
the two-form source J^ x , causes an electromagnetic field = \e i i VK ,\F K \ that 
is now Meissner screened over a length scale A. The case p, = t, together with 
the definition B n = \e na bF a b reduces to Eq. (jT7|) . 

Summarizing, we have shown here that the action Eq. (f2l?)) can be postu- 
lated, from which the correct equations of motion as introduced in Section [3] 
directly follow, without ever mentioning the gauge potential A^. One trades 
in the Bianchi identities Eq. (|19l) for a set of gauge transformations Eq. ((32)) . 
This action is only meaningful in the absence of monopole sources, but is very 
appropriate when considering two-form sources such as Abrikosov vortices. In 
the case that the penetration depth A becomes infinitely large, the field strength 
F^y recovers its status as a gauge field. This would correspond to the Coulomb 
phase of two- form sources, as opposed to the Higgs phase that is always realized 
in superconductors. 

As a final note it should be stressed, that although the vortex source is 
intrinsically dipolar in nature, the equations stated above are not generally valid 
for any dipole source. Here, the direction of the vortex line is always parallel 
to the dipole moment. If one should instead consider for instance a string of 
ferromagnetic material with moments not along the string, one must revert to 
the omnipotent regular Maxwell equations. 

For the reader familiar with differential forms, we have included an appendix 
repeating these considerations in metric-independent language, valid in any spa- 
tial dimension higher than 2. 



5 Vortex duality in charged superfluids. 

We shall now rigorously derive the coupling of Abrikosov vortex sources to the 
electromagnetic fields, starting from the action describing a superconductor in 
3+1 dimensions. This follows the same pattern as the first steps of the well 
known vortex (or Abclian-Higgs) duality in 2+1 dimensions, where by Legcndrc 
transformation the vortex-vortex interactions are expressed in an effective elec- 
trodynamics sourced by the vortex sources. In a recent publication we demon- 
strated how to extend this to 3+1 dimensions [l6|, in a language that is very 
closely related to the present context, also resting on two- form gauge theory. 
However, in this earlier work we concentrated on the neutral superfluid and here 
we will go one step further by coupling in electromagnctism and integrating out 
the condensate gauge fields, ending up with an effective action describing the 
electromagnctism of vortices. 

Our starting point is the partition function Eq. (|13[) . To keep the equations 
readable, we will transform to dimensionless units denote by a prime (which we 
suppress when matters are unambiguous), 

, 1 , 1 , c , ae* , ha 2 1 h 

& = T^i x m = ~ x m, t = —t, A = ——A^, p = —Ps, —. = -y- 

h a a p n m*c pi Poce* 

(36) 

Here a is a length scale relevant in the system, for instance the lattice constant. 
We will assume the relativistic limit c p h = c; later we shall return to dimen- 
sionful quantities and it will become clear that the phase velocity is playing an 
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essential role for the description of the non-relativistic vortices. The partition 
function in these dimensionless units reads, 

ifd 4 x C 



Z = J VpVA^(A fl )e 1 J ax L , (37) 
C = -^% > -\p(d ll <p-A ll f. (38) 

Now we perform the dualization procedure. A Hubbard-Stratonovich trans- 
formation of Eq. (|37)) leads to, 

Z = J Vw^VpVA^(A^e^ c ^\ (39) 

^dual = -^flv + ~ ~ (40) 

Here w M is the auxiliary variable in the transformation, but it is actually the 
canonical momentum related to the velocity d^ip, which can be found as, 

dC 

w ^ l = ~d(d^) =p ( dflip ~ A ^ 1 (41) 

and is related to the supercurrent as = %- . If one integrates out the field 
from Eq. (|39[) . one retrieves Eq. (|37|l . In the presence of Abrikosov vortices, 
the superconductor phase ip is no longer everywhere single-valued. Therefore 
it is separated into smooth and multi- valued parts ip = (^smooth + Pmv- The 
smooth part can be partially integrated yielding, 

Z = Jvw^ Vp smooth Vpuv VA^ F{A„) e 1 / , (42) 

£dual = --7~ F lv + -W^Wf, + ^smoothd/^ - W^^MV + W^A^. (43) 
4p p 

The smooth part can now be integrated out as a Lagrange multiplier turning into 
the constraint d^w 11 = 0, the supercurrent continuity equation. This constraint 
can be explicitly enforced by expressing as the curl of a gauge field. In 3+1 
dimensions, this gauge field is a two-form field, 

vf = e^ KX d v b KX . (44) 

We can now substitute this expression in the partition function; the integral 
over the fields u> p is replaced by one over b KX , as long as we apply a gauge fixing 
term T(b K \) to take care of the redundant degrees of freedom. Since the gauge 
field is smooth it can be partially integrated to give, 

Z= fvcp M y'DA^^(A^Vb KX ^(b KX )e i S c ^, (45) 
£duai - -7^1 + -{^ VKX d v b RX f - b KX e KX »»d v d^MV + b^e^d^. 

4:fl p 

(46) 

Here we recognize the definition Eq. (jllj) of the vortex source, 

JJa = e^^^lfttv, (47) 
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and we have derived the dual partition function, 



Z = Jvj^VA^TiA^Vb^Tib^e^^, (48) 

£duai - ~Flu + -(^ KX d„b KX ) 2 - b KX J^ x + b^e^dvAp. (49) 
4/1 p 

The interpretation is as follows. The vortex sources interact through the 
exchange of dual gauge particles b K \ coding for the long range vortex-vortex 
interactions mediated by the condensate. The gauge field b K \ couples as well to 
the electromagnetic field A p . Integrating out the electromagnetic field will lead 
to a Meissner/Higgs term ~ showing that the interaction between vortices is 
actually short-ranged in the superconductor. However, we are instead interested 
in how the electromagnetic field couples to the vortices themselves. Therefore, 
we shall integrate out the dual gauge field b K \. 

The first step is to complete the square in b K \. The kinetic term for b K \ is 
proportional to, 

- b KX ^ v d v e paa »d a bP° = -b KX (S^8 2 - d K d^b^ = -b KX Qo 1K %X- (50) 

Here Gq 1 ^^ is the inverse propagator. However, this expression cannot be in- 
verted (the same problem arises in the quantization of the photon field). We 
can solve this by imposing the Lorenz gauge condition d K b K \ = 0. Then the 
inverse propagator is simply Qq 1K = 5 Kf *d 2 , and its inverse is Qo Kll = <5 K/i ^2- 
Now we can complete the square, 

£duai = \{b^ £ J K V A + e KXufl d v A») ( - -) {b^ - £ J VKA + e^d p A a ) 

~\{-J^ + W™«) ( - £) ( - J VKA + e^d v A») - ±-F%,. 

(51) 

Then we shift the field b K \ — > b K \ + -§?J^\ — ^nXv^d 11 A^ and integrate it out 
in the path integral to leave an unimportant constant factor. Expanding the 
remaining terms leads to, 

C A , - I T V P_ jV kX , 1 cf, ah J^ KXpa * A _ jV K \v^ A _J_p2 

i-dual — ^ J kA q 2 J +2 ^"^ Q2 °P A <y P J K \ e d 2 iu ^ 

= l^ VKX - \p^A, - P Jl^%A, - (52) 

In going to the second line we have performed partial integration on the second 
term and invoked the Lorenz gauge condition d^A^ = 0. We can immediately 
read off the physics encoded in this action: the first term describes the core 
energy of the vortices and we shall not need it in this work; the second term is 
the Higgs mass (including Meissner) for the electromagnetic field; the third term 
is the coupling term between the electromagnetic field and the vortex source. 
This last term looks rather awkward given the derivatives in the denominator. 
This could signal that the coupling is non-local but that is not the case here. 
The origin of this coupling follows from the notions presented in section 21 it is 
not the gauge potential A p but rather the field strength F pv itself that couples to 
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the vortex source. We can confirm this expectation by computing the equations 
of motion, 

-d„F^ + pe^ x ^ jv - p A v = 0. (53) 
Acting with Cap-yud 1 on this equation, one obtains, 

— e aMv dW^ + e aMv e»™ x ^jV x - e ahv d^ A v = (54) 
pp u 

Using F M „ = d^A v — d^A^ one can see that from the first term only e a fs t _ lu d 2 F fJ - 1 ' 
survives. Also, using d K J^ x = one can see that ^ a ^vf liUKX d 1 d il J^. x = d 2 J^p, 
cancelling the derivatives in the denominator. Altogether we find, 

^e^a 2 F^ - \e aP , u F^ = -jv . (55) 

This is the same result as Eq. p5[). As we announced earlier, we have derived 
here with a completely controlled procedure the dimensionless version of Eq. 
(jTTJ) , describing the interactions between the vortices and electromagnetic fields 
inside a relativistic superconductor. Departing from this result we will derive 
in the next section various physical consequences. Summarizing this section, 
by dualizing the Ginzburg-Landau action for the superconductor, Eq. (|37p got 
reformulated in terms of the vortex currents Eq. (|47|) as the active degrees of 
freedom, that interact via the effective gauge fields parametrizing the rigidity of 
the superconductor. The latter were integrated out to obtain the direct coupling 
of the vortices to the electromagnetic field, leading eventually to the concise 
equations of motion Eq. (f55"j) . Although this strategy is well known dealing 
with vortex 'particles' in 2+1 dimensions we are not aware that it was ever 
explored in the context of the electrodynamics of vortices in 3+1D. Surely, the 
derivation presented in the above is in regard with its rigour and completeness 
strongly contrasting with the rather ad hoc way that the problem is addressed 



in the standard textbooks [15J, eq.(5.13) 



6 Vortex electrodynamics. 

In order to establish contact with the physics in the laboratory all that remains 
to be done is to break the Lorentz invariance, doing justice to the fact that 
the phase velocity of the superconductor as introduced in the first paragraphs 
of Section [3] is of order of the Fermi velocity of the metal and thereby a tiny 
fraction of the speed of light. Subsequently we will analyze what the physical 
ramifications are of our "Maxwell equations for vortices" . 
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The non-relativistic version of the vortex action Eq. (|52")l is, 



r-——n fV 1 jV 11 jV °ph tV 

*2 *2 

e ,.n e 



2m*Cp h 2?n* 

1 .32 T a2 [ J^b € abtm(dtA m + 8 m V) + — J^etamn<9 m j4 n ] 

ph 

+ ^-jCfleA. + 5„^) 2 - -^-(d m A n - d n A m f. (56) 

Varying with respect to A v and acting with e Q ^ 7l/ (9 7 and imposing current con- 
servation d K J^ x = will lead to the correct equations of motion. However, the 
easiest way to obtain the non-relativistic versions of the equations of motion 
Eq. ([TB")) is to vary Eq. (TlH) directly with respect to V and A n respectively, 

-%A 2 d„£„-V = 4^, (57) 
cr e* 

— A 2 ^-3 t £' n + X 2 e nmk d m B k + A„ = 4 9 «¥>- ( 58 ) 
cr e* 



Here A = ^/ ^*i p is the London penetration depth. Now we operate on 



the first equation by d m = \e m tab^abrtdr , and on the second by S mn d, 



9 ^traab tabtndt and t tamn d m respectively to obtain, 



2 



c /i 1 

5- A d m d n E n — 9 m V = — c p h-e mta b J ab , (59) 

e* 2 

— A 2 ^- d^E m — A 2 <9 2 i9 t A m + 9 t A m = — c p h-e fma b J^, (60) 
A2( V 2 _la 2 )B Q -B a = -AjV (61) 

For the last equation we used the Maxwell equations V x E = — <9 t B and V • B = 
0. It is equal to the one we found before in Eq. (fTT)) . obviously, since there the 
temporal terms do not play a role. 

For the equations for the electric field is it useful to choose the Coulomb 
gauge V • A = 0, and separate the electric field in longitudinal and transversal 
parts: E = E L + E T , where V x E L = and V • E T = 0. In the Coulomb gauge 
we see from the definition E = — W ~ d t A that E L = -VV" and E T = -d t A. 
We can subtract the first equation above from the second to obtain, 

A 2 ( ^d^E m + V 2 i?^ + -^-V 2 E^) — E m = —Cp^etmabJab- (62) 

c c e 

Hence, as in the case of the Maxwell theory for non-relativistic matter one 
finds instead of the highly symmetric relativistic result Eq. (|16p two equations 
of motion that are representing the spatial- (magnetic) and temporal (electrical) 
sides of the physics, Eq. (jBTj) and Eq. (j62|) . One notices that the first 'magnetic' 
equation is quite like the relativistic one while the 'electrical' equation is now 
more complicated for reasons that will become clear in a moment. 
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The factor c p h on the right-hand side of the electric equation is due to our 
convention of rescaling the time derivative to having units of 1 /length in the 
definition of J^ x . Thus all components of J^ A have dimensions of a surface 
density, and multiplying by a velocity is necessary to end up with a current 
density. The sign difference on the right-hand side between the electric and 
magnetic equations is related to the continuity equation ^i_9 t = — d m J^ b . 

To grasp the content of these equations, one should compare the magnetic 
equation Eq. (|6ip with the standard form[15j, eq.(5.13)], 

\ 2 V 2 B a -B a = -$ <5i 2) (r), (63) 

Here $o = 2nh/e* is the flux quantum. The factor of 2n is associated with the 
definition of J v as in Eq. (|11[) . Our treatment automatically takes dynamics 
into account in the form of temporal derivatives. Otherwise, the correspondence 
is complete. We have indeed exactly recovered the well-established vortex equa- 
tion of motion. 

The equation for the electric field (1621) looks more involved, but this can be 
made more insightful by writing the equations for the longitudinal and transver- 
sal parts separately, 

^((^V 2 - M 2 )^ - El = ^c ph el iab Jl, (64) 
A 2 ((V 2 - ^d 2 t )El - El = ^c ph ej mab jv. (65) 

The labels on the e-symbol denote that they include a longitudinal or transversal 
projection. 

We can now read off the following physical relations: 

1. Meissner screening: from Eq. (|61[) in the static limit d t — > 0, a vortex 
line sources a magnetic field, that falls off in the superconductor with a length 
scale A, the familiar Meissner effect. 

2. Thomas-Fermi screening: from Eq. (|64p one infers that the longitudinal 
(electrostatic) electric field penetrates up to a much smaller length ^-A, which 
is the Thomas-Fermi length (c « 300c p h). This just amounts to the well known 
fact that the electrical screening is the same in the metal as in the superconduc- 
tor. Notice that this length scale is obtained without referral to the electrons 
in the normal metal state as in the textbook derivation. 

3. Dynamic Meissner screening or the Higgs mass: taking into account the 
time-dependence, Eq. (|6"Tj) and Eq. (|65[) show that the transversal photon parts 
of the fields are screened not only in space, but also in time with characteric 
time scale — . This is just the familiar statement that the two propagating 
photon polarizations in 3+1 dimensions acquire a "Higgs mass" ~ inside the 
superconductor . 

4. Electrical field of a moving vortex and the Nernst effect: disregarding the 
dynamical term in Eq. (|62[) . one is left with 

E m = ——Cph^tmklJkl- (66) 

Recall from section [2] that we had interpreted as the flow or velocity in the 
^-direction of a vortex line in the / direction. Since we know that one vortex 
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(a) A vortex in a Joscphson junction between two (b) Vortex. Geometry of the electric field E 

superconductors (grey); it has no normal core, generated by a vortex line parallel to the mag- 

The magnetic field B is along the vortex; any elec- netic field B and moving with a speed v. This 

trie field across the junction causes the vortex to phenomenon related to the Lorentz force follows 

move in the perpendicular direction. Such mo- directly from the vortex equations of motion, 
tion induces electromagnetic radiation that may 
escape to the outside world. 

Figure 3: Additional vortex configurations. 



line carries a magnetic flux of $o = 27Tp-, we can write p-c pn J^i = v iB®, where 
B° denotes the field associated with one quantum of flux, and vi = c pn ez is the 
velocity. In practice there is always a drag force that greatly slows down the 



vortices. Still, Joscphson vortices that do not have a normal core (Fig. 3(a) ) 
may achieve this large speed. With this interpretation, (pj^j) reads, 

E = -v x B°, (67) 

which is precisely the known result [l7| for the electric field generated by a 
vortex moving in a magnetic field B° (Fig. |3(b)[ ). When the motion is caused 
by a temperature gradient this is responsible for the large Nernst effect of the 
vortex fluid. 

5. AC Joscphson relation: another interpretation of Eq. is found by 
inserting the definition of the vortex current, J^ b = e a bt n -^dtd n ip, taking m as 
the longitudinal direction and neglecting the higher derivative terms. In this 
case, 

d m V= ^d t (d m ip). (68) 

e* 

Here the left-hand side is the potential difference, and the right-hand side is the 
time derivative of the superconducting phase difference. This is exactly the AC 
Joscphson relation. The full equations Eq. (|6"!2)) reveal also that the induced 
electric field is screened inside the superconductor. 

6. Moving vortices as radiation sources: in the same spirit, the moving vor- 
tex is also inducing dynamic transversal fields according to Eq. (|65p . In other 
words: moving vortices radiate Q. But since the field is Mcissner screened, it 
is very hard to detect this radiation. All our results also apply to Josephson 
vortices (line vortex solutions in a Josephson junction between two supercon- 



ductors parallel to the interface, Fig. 3(a) I, which differ only in the regard that 



they do not have a normal core. There is much recent interest in radation from 
(arrays of) Josephson junctions, see e.g. [l8|. Since inside the junction the field 
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is not expelled by Meissner and metallic screening, the radiation may escape to 
the outside world. In this literature one finds the following result[19j, eq.(13)], 



- A 2 V 2 A + A = — V0. 



e* 



(69) 



Here A differs from A because of a special geometry. Compare this with a result 
that follows from Eq. ((65]) . 



confirming Eq. (|69[) but showing in addition how to take care of a possible time 
dependence of the photon field. 

Summarizing, to the best of our knowledge we have addressed all known 
electrodynamical properties of vortex matter departing from a single action 
principle. 

7 Conclusions. 

We are of the opinion that our action principle for vortex electrodynamics Eq. 

rcsp. (fSB")) and the associated "vortex-Maxwell" equations Eq. (|6"Tj) 
and (f62|) deserve a place in the textbooks on the subject. In contrast with 
the clever but improvising discussions one usually finds, our formulation has 
the same 'mechanical' quality as for instance the Landau-Lifshitz treatise of 
electromagnetism. One just departs from the fundamentals, to expose the con- 
sequences by unambiguous and straightforward algebraic manipulations that are 
worshipped by any student of physics. A potential hurdle is that one has to get 
familiar with the two-form gauge field formalism, but then again this belongs 
to the kindergarten of differential geometry and string theory. 

Our analysis also reveals the origin of the peculiar nature of this vortex 
electrodynamics. The realization that it is in fact governed by a two-form gauge 
structure amounts to an entertaining excursion in the fundamentals of gauge 
theory itself, nota bene associated with the superficially rather mundane and 
technology-oriented vortex physics, at least when viewed from the perspective of 
fundamental physics. In fact, our pursuit started in the quite different subject 
of the dual description of Bose-Mott insulators in terms of superconducting 
vortex matter where quite similar issues arise. In this context one does not get 
anywhere with physical intuition and one just needs full mathematical control 
to make any progress. Therefore we decided to inspect these matters first in the 
more familiar present context. 

On the practical side, as we implicitly emphasized in the last section our ap- 
proach offers a unified description of the electrodynamics of vortices. Although 
we got as far as recovering the known physical effects in terms of special lim- 
its of our equations, we sense that there is a potential to use them to identify 
hitherto unknown effects and perhaps to arrive at a more complete description 
of the electrodynamics vortex matter. Being well aware of the large body of 
knowledge of this large field in physics, we leave this as an open question to the 
real experts. 



d t [-A 2 ((V 2 



M 2 )A + A]=9 f [-§V^], 



(70) 
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name field ?-form 2+ld 3+ld representative 



in d=3+l 



electric field 


E 




1 


1 


1 


E x dx 


dielectric current 


D = 


= e * s E 


d-2 


1 


2 


D x dy A dz 


magnetic field 


B 




2 


2 


2 


B x dy A dz 


magnetic intensity 


H = 


= (i* s B 


d-3 





1 


H x dx 


charge density 


P 




drl 


2 


3 


p dx A dy A dz 


current density 


J 




d-2 


1 


2 


J x dy A dz 


covariant current 


j = 


p + J A dt 


drl 


2 


3 


j x dy A dz A dt 


field strength 


F = 


= B + E A dt 


2 


2 


2 


F xy dx A dy 


gauge potential 


A 




1 


1 


1 


A x dx 


vortex source 


J v 




d-2 


1 


2 


Jxy A dy 


Lagrangian density 


£ 




d 


3 


4 


C dtAdxAdyAdz 



Table 1: Electrodynamical quantities in differential forms. Here e and p are the 
electric permittivity and the magnetic permeability, and * s is the spatial Hodge 
dual. Other factors of c are suppressed. Minus signs are subject to convention. 
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A Electrodynamics with differential forms. 

For the reader familiar with the mathematical language of differential forms, we 
present the electrodynamics of vortex sources for any dimension d = D+l higher 
than 2. For our purposes, a differential form can be thought of as something that 
can be integrated over; in other words: it is a density function combined with the 
integrand. For instance, the electric field is a 1-form E = Eidxi = E x dx+E y dy+ 
E z dz. Higher forms are always obtained through the wedge product a Ab, which 
is the antisymmetrization of the tensor product of a and b. Another common 
operation is the Hodge dual *a of a, which turns an n-form into a (d— n)-form. 
For instance in three spatial dimensions *E = E x dyAdz+E y dzAdx + E z dxAdy. 
For a pedagogical introduction to differential forms in Maxwell electrodynamics 
see [20l | . 

In the familiar case of d = 3 + 1, a 1-form is a line density or "field inten- 
sity" like the electric field; a 2-form is a surface density or flux density like the 
magnetic field; a 3-form is a volume density like the charge density. Confusion 
may arise when it is not immediately clear whether an object is an n-form or 
a d — ?i-form, which is important for generalization to other dimensions. We 
distinguish the regular Hodge dual * from the spatial Hodge dual * s , where the 
latter does not involve the temporal dimension. The exterior derivative opera- 
tor is d = d< A + dxiA, and the one with only spatial components is 

ds = J2i~£r d^A. The Leibniz rule is d(a A b) = da A b + (— l) r a A db, where a 
is an r-form. This can be used for partial integration. 
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In tableQ]we have listed the differential forms of the relevant fields. Some of 
these definitions seem perhaps unfamiliar. In particular, we are used to thinking 
of the magnetic field as a vector field; however, its solenoidal nature is typical 
of a two-form. This becomes even more clear when it is expressed as the curl of 
the vector potential B = d s A, which holds in 3 + dimensions. Also the current 
density J is naturally a flux or a 2- form, but its generalization is as a d — 2-form. 
One way to see that this must be so, is to write down the continuity equation 
in differential forms, 

<9 t p + V-J = -> (d t p + d s J) A dt = dj = 0. (71) 

The current density appearing as a vector field for instance in Ohm's law, J = 
crE is actually the spatial Hodge dual of J. 

We shall now write down the familiar expressions of Maxwell electrodynam- 
ics. The Lagrangian density is a spacctimc volume density. All terms must 
therefore combine into d-forms. The field strength is F = dA. From this defini- 
tion it is clear that the gauge transformations A^A + d£, with £ any 0-form, 
leave the field strength unchanged, since d 2 = 0. The field strength is contracted 
with its dual to obtain a <i-form in the Lagrangian. The sources couple to the 
gauge potential (this is another reason why the source is a d — 1 form) . The 
Maxwell action is then, 

S = y — F A *F + A A j. (72) 

The second term is also invariant under the same gauge transformations, pro- 
vided that dj = 0, the continuity equation. The Euler-Lagrange equations are, 

resulting in the inhomogeneous Maxwell equations, 

d * dA = d * F = -j, (dpF/xv = —Jv)- (74) 

Applying the exterior derivative on this equation directly leads to the continuity 
equation, since d 2 = 0. Similarly, from the definition F = dA it immediately 
follows that, 

dF = 0, (75) 

which are the homogeneous Maxwell equations, or in this context rather the 
Bianchi identities. 

Now let us repeat the reasoning of section @] In the absence of monopolc 
sources J, we have both d * F = and dF = 0. This implies that the field 
strength has become "pure gauge" . The first of these equations still holds when 
we add any 1-form £ as *F— ^ *F + d£. The original Bianchi identities are not 
invariant under these transformations. The dual field strength *F turns into 
a gauge potential, and is accompanied by its own field strength K = d * F, 
which contracts with its dual in the Lagrangian. The field strength can couple 
to a d — 2-form source, which we anticipatingly denote by J v , provided that 
d * J v = 0. Indeed, 

F A J v -> F A J v + *d£ A J v = F A J v - £ A d * J v = F A J v . (76) 
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The second step is achieved by partial integration, and the last equality holds 
if the vortex current is conserved, d * J v = 0. The action for vortices directly 
sourcing the field tensor is, (with G = *F), 

S= /"-KA*K + FAJ V = y*-KA*K + GA*J v (77) 

Variation with respect to G leads to, 

*d*dG = -J v . (78) 

This equation corresponds to e K \ li vd 2 F fll , = —JY\ as in eq. (|55p . but is valid in 
any dimension. The addition of a Mcissner term results in 



S = J — K A *K — G A *G + G A *J . (79) 

and, 

*d*dG-G = -J v . (80) 

This is the equation of motion for d — 1-dimensional superconductors, which 
have d — 2-dimensional vortex worldbranes J v . 
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